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I. INTRODUCTION 

The Dicke model describes the dynamics of A" identi- 
cal two-level atoms interacting with a quantized three- 
dimensional electromagnetic (EM) field Under cer- 
tain conditions the model predicts that the atoms inter- 
act with the quantized EM field collectively, giving rise 
to the widely studied phenomena of superradiance and 
subradiance 0, H[ . In free space ideal superradiance and 
subradiance take place in the so called small sample limit, 
i.e., when the atoms are so close to each other that one 
can ignore any effect resulting from their different spa- 
tial positions. In this case the atoms are indistinguishable 
with respect to their emission and absorption properties; 
hence, the presence of equivalent paths through which the 
emission process may occur gives rise to fully construc- 
tive (superradiance) or destructive (subradiance) inter- 
ference. 

Ideal superradiance or subradiance in free space is very 
difficult to observe in the experiments since it requires 
that the atoms are placed in a regular pattern within a 
sample smaller than the wavelength of the EM field they 
interact with (small sample case). The requirement of 
a regular pattern is due to the presence of the dipolc- 
dipole forces that would otherwise break the symmetry 
under permutation of any two atoms necessary to observe 
supcrradiant-subradiant behavior. Such a regularity can 
be achieved, e.g., with trapped- ion crystals [4[ or atoms 
in optical lattices ||. In these systems, however, the sep- 
aration between the particles is typically larger or on the 
same order of magnitude than the resonant wavelength 
(large sample case). In the large sample case, cooper- 
ative effects still occur but the subradiant state is not 
completely decoupled from the dynamics. Indeed, par- 
tial subradiance and superradiance have been observed 
with trapped ions Q. 

A way for relaxing the requirement for configuration 
regularity is to place the small sample in a cavity res- 
onator. In this case, indeed, due to the Purcell effect, 
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the cooperative atomic behavior can be observed at much 
lower atomic density than in free space, making the van 
der Waals dephasing caused by the irregular atomic con- 
figuration negligible Q. Experiments observing super- 
radiance in the small sample case in a cavity have been 
performed with Rydberg atoms Q, giving results in a 
very good agreement with the predictions of the single- 
mode superradiance theory. In this experiment, all of the 
atoms are equivalently coupled to the quantized mode of 
the EM field (homogeneous case). 

Recent advances in ion-cavity QED experiments make 
it possible to confine arrays of ions inside an optical cav- 
ity in a regime in which the width of their wave packet in 
position space is smaller than the wavelength of the cav- 
ity mode they interact with (Lamb-Dicke regime) [1, [(| . 
Moreover, it is possible to accurately manipulate the po- 
sition of the single ions with respect to the intensity pro- 
file of the standing cavity mode, thereby allowing us to 
change the strength of the coupling between each ion and 
the quantized EM field. 

It has been demonstrated theoretically that, when the 
atoms are coupled with different strengths to the EM 
field, ideal superradiance or subradiance can still occur, 
depen ding on the particular spatial distribution of the 
atoms [HI [III [IH ■ However, no experiments have up to 
now confirmed these predictions by the inhomogeneous 
Dicke model. Very recently, an important step in this 
direction has been achieved at the University of Aarhus, 
where a collective strong coupling between an ion crys- 
tal and a cavity mode was observed In this paper, 
we investigate in detail how the inhomogeneous single- 
mode Dicke model (or Tavis-Cummings model [l4j ) can 
be realized in the ion-cavity QED context and the con- 
ditions under which subradiance and superradiance can 
be observed. 

Besides the importance in the study of fundamentals 
of quantum theory, the realization of the Dicke model 
and the generation of the subradiant state play a crucial 
role in quantum information technology and quantum 
communication. Indeed, arrays of ions are ideal candi- 
dates for quantum registers and their controlled interac- 
tion with photons allows us to realize atom-light quantum 
interfaces [l5[ and to distribute entanglement to different 
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nodes of quantum networks. The importance of the sub- 
radiant states in this context stems from the fact that 
they arc robust entangled atomic states since they are 
completely decoupled from the EM field. 

The aim of this work is to discuss a realistic setup that 
is able to show the collective behavior of trapped ions 
in a cavity. In particular, since in the experiments per- 
formed so far the ions are coupled to the EM mode via 
a Raman scheme in a A-configuration, we will include 
the entire level structure, which is important in order to 
understand the decohering role of the spontaneous emis- 
sion from the upper and essentially unpopulated level. 
We will also include cavity losses in order to study in de- 
tail the deviation from the ideal cooperative Dicke model 
and to identify the parameter regions in which such de- 
viations are as small as possible. 

In fact, during the last two decades, several theoretical 
papers have discussed issues such as entanglement gener- 
ation, preparation of nonclassical states, or realization of 
quantum gates in the ion-cavity QED context assuming 
that the conditions to realize an ideal Tavis-Cummings 
model were met [H M, M M M M M M M M, M- 
Thus, cither the spontaneous emission or the cavity losses 
(or both processes) are usually neglected [l?], O US HI] • 
Concerning spontaneous emission, for example, the as- 
sumption is made that the emission rate is much smaller 
than the cavity coupling constant 0, H3, HH, [H, HH- 
However, this condition is not met in the ion-cavity QED 
experiments @, Q. Furthermore, as we will demonstrate 
in this paper, if one deals with simplified atomic level 
structures [3, [H, HE H3] , it is not possible to single out 
those regions in parameter space for which the systems 
of trapped ions behave collectively. 

In this paper, we will take both the cavity losses and 
the spontaneous emissions into account and employ A- 
type schemes to describe the ions and to identify the ex- 
perimental conditions under which the coherent dynam- 
ics predicted by the single-mode Dickc model is dominant 
with respect to losses and decoherence. This will also 
allow us to present realistic protocols for entanglement 
generation and to discuss ways to optimize the generated 
entanglement using specific features of the trapped-ion 
system, such as the ability to manipulate in a controlled 
way the relative coupling between the ions and the cavity 
field. 

The structure of the paper is the following. In Sec. |TT] 
we review the properties of the inhomogeneous single- 
mode Dicke model. In Sec. IIIII we present the Hamilto- 
nian for two ions in a cavity and we make the connection 
to the Dicke model by deriving an effective model de- 
scribing the dynamics under realistic experimental con- 
ditions. Section IIVI is devoted to the description of the 
experimental proposal to observe subradiance and verify 
the inhomogeneous Dicke model. Furthermore, in Sec. [V] 
we explore another way to optimize the entanglement 
generation by using off-resonant transitions. Finally, a 
summary of the results and the conclusions are given in 
Sec. EH 



II. INHOMOGENEOUS SINGLE-MODE DICKE 
MODEL 

A. Ideal cavity 

The single-mode Dicke model, or Tavis-Cummings 
model, is the simplest quantum-mechanical model de- 
scribing collective effects such as supcrradiancc and sub- 
radiance in cavity. It describes the quasi-resonant inter- 
action between N identical two-level atoms and a single 
quantized cavity mode. The Tavis-Cummings Hamilto- 
nian is 

N 

+ (a^aV^ + a«*a<4>')) , (i) 

where uc and oja arc the frequencies of the cavity mode 
and the atomic transition, respectively; a and are 
the annihilation and the creation operators for the cav- 
ity mode; and a U) = |0^)(1^| and = (oi j) ) f are 
the lowering and the raising operators for the jth atom, 
|0w) and |D-") being its ground and excited states, re- 
spectively. Finally, a® is the coupling strength of the 
jth atom with the cavity field. Inhomogcncity of the cou- 
pling strengths originates from different relative positions 
of the atoms with respect to the intensity profile of the 
standing EM mode supported by the cavity resonator. 

This model assumes that the cavity is ideal, as pho- 
ton escape is not taken into account, and that atomic 
spontaneous emission from the excited to the ground 
state is negligible. The model also neglects the atomic 
motion as well as recoil effects due to the absorption 
and subsequent re-emission of a photon by the atoms. 
Moreover, the dipolar coupling of the atoms and the EM 
field is expressed within a rotating wave approximation 
(RWA), thereby suppressing the non-energy-conserving 
terms. Finally, it implicitly assumes that the coupling 
between the atoms and the cavity mode does not change, 
i.e., that the atoms are kept at fixed positions. While the 
RWA has been proven to work extremely well in optical 
experiments, all other assumptions need further consid- 
eration. In the following sections we will examine them 
in detail for the ion-cavity QED setup. 

Using a suitable canonical transformation it has been 
shown that, when only one excitation is present in the 
total system, the N atoms interacting with the quan- 
tized field mode according to Eq. (TT|) cooperate in such 
a way that only one collective atomic mode (superradi- 
ant state) is coupled to the field [llj. Consequently, the 
energy exchange between the atoms and the field can be 
completely suppressed if the only field-coupled collective 
mode is unexcited. 

For simplicity, wc will from now on focus on the N = 2 
case sketched in Fig. [TJ and we will denote the energy 
cigenstatcs for the free ions as ja^M 2 ') = \a^) g) \b^) 
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FIG. 1: (Color online) Two binary quantum objects interact- 
ing through a quantized electromagnetic mode supported by 
a cavity resonator. The dynamics of such an ideal system is 
described by the Dicke model. 



(with a, b = 0, 1) and the corresponding Fock states of 
the cavity mode as \n^), where n = 0, 1, . . . . The time 
evolution generated by Hjj is easily obtained explicitly. 
For a cavity initially prepared in the vacuum state, and 
in the presence of only one atomic excitation, the time 
evolution of the amplitudes cio(i) and Coi(i) to find the 
ions in the states and \0^1^), respectively, is 

given by 



cio(*) 
coi (t) 



\r^\ 2 + \r^\ 2 £(t)]c w (0) 
. r (DV 2 )[l-£(t)]c 01 (0), (2) 

. r (iy 2 )*[i-£(t)] Clo (o) 

|rW| 2 + |rW| a f(t)lcoi(0). (3) 



In the equations above the relative coupling strengths 



arc defined as 



Si) , 



lorl, where \olt\ 



y/\aW | 2 + |a( 2 ) | 2 is the total coupling strength, and 

£ (t) = e iSt l 2 r 
where 5 



cos l— J-iT sm l— ) 



(4) 



u>c is the detuning and Q v = 
- 5 2 is the vacuum Rabi frequency. Note that 
A 2 ^ are not independent parameters, since 



n 1 ) and 

| r (2)| = | r (l)|2_ 

The subradiant 
are 



B. Non-ideal cavity 

We now proceed to generalize Eq. (QJ to the case of 
a lossy cavity. The imperfect reflectivity of the cav- 
ity mirrors and consequent leakage of photons causes a 
Lorentzian broadening of the spectral line correspond- 
ing to the mode supported by the ideal cavity. Accord- 
ingly, the microscopic atom-field interaction should now 
take into account a continuum of modes described by a 
Lorentzian distribution peaked at the central cavity fre- 
quency u>c ■ For the sake of simplicity, and in view of the 
discussion in the ion-cavity QED context, we restrict our 
attention to a one-dimensional cavity model. Namely, we 
neglect the coupling with all the EM modes other than 
the ones supported by the lossy cavity. In the rotating 
wave approximation, the Hamiltonian is given by 



H 



i 



JV 



k a k 



u>act + <J_ 



3 = 1 



N 



]T [to sin (^r x(j) ) a W- + H - c - 



(7) 



k 3 = 1 



where a k and a k are the annihilation and the creation 
operators of cavity photons of frequency u> k , respectively. 
Above, we have assumed that all the atoms have the same 
electric dipole moment, which has been incorporated in 
the coupling constants g k , and we indicate with x^' the 
position of the atoms along the cavity axis. In the fol- 
lowing we will assume that each atom is kept at a fixed 
position inside the cavity and that they are all well lo- 
calized, i.e., the spread of their wave function in position 
space is smaller than the wavelength of the central cav- 
ity field mode: Ax^ -C c/ujc- Since all the significantly 
contributing modes are close to the central mode (of fre- 
quency ujc), we have 



sin — x yJI 
V c 



sm 



(8) 



and the superradiant states and Eq. takes the form 



|^_>= r (2)|i(i) ( 2 )>-r«|0«l( 2 )>, 



(5) 
(6) 



1 N 



3 = 1 



and, in this case, they are position dependent through the 
relative coupling strength parameters A x > and A 2 \ As 
one can see directly from Eq. ([T]), the state |Cr c ') is 

an eigenstate of the Tavis-Cummings Hamiltonian with 
eigenvalue \<^c + ^>A- Therefore, when the atoms are 
prepared in this state, they arc completely decoupled 
from the cavity field and the system does not evolve 
at all. In the case of equally strong couplings, i.e., for 
| r (i)| = | r ( 2 )| = 1/V2, the subradiant and the superra- 
diant states coincide with the maximally entangled Bell 
states. In general, however, these states are not maxi- 
mally entangled. 



N 

J2[ X ^J2g k al + H, 



(9) 



3 = 1 



with x = i sux(locx^' / c). In the continuum limit the 
sum over the fc-modes is replaced with an integral 



E 



9k\ 



dwJ(uj), 



where J(oj) is the reservoir spectral density. As men- 
tioned above, we assume a Lorentzian distribution for 
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the spectrum of the field inside the cavity; therefore, we 
take a spectral density of the form 



J(w) 



W 2 



2tt 



(w - ujcY + («/2) 5 



(10) 



where the distribution is characterized by its full width at 
half maximum value n and by a normalization parameter 
W 2 = Jdio J(ui). Hence, K describes the cavity losses and 
W describes the total coupling strength. 

We focus again on the two-atom case, i.e., N = 2, and 
we consider the situation in which only one excitation 
is present in the total atoms-field system. Starting from 
the Hamiltonian ([7]) and using the Lorentzian spectral 
density (JTUJ) , it is possible to derive an effective master 
equation 

— ^ = — i [Hd, q] — — [aftag + gefta — 2aga l ] (11) 

for the dynamics of the atoms and the cavity mode of 
frequency u>c [28| . Here, a and aft are the annihilation 
and the creation operators for the central cavity mode, 
which is damped at rate k, and the coherent dynamics 
is generated by Hp in Eq. (fT]), where the coupling con- 
stants are identified as a* = x W- From the exact 
solution of the effective master equation (fTTj) , one can 
obtain the state of the atomic system by tracing out the 
cavity degree of freedom: p(t) = trc[g(i)]- 

After performing the trace, and for an initially empty 
cavity, the problem can be solved exactly. In the 
ordered basis \1^0^), |0«l( 2 )}, loW*}}, 

the atomic density matrix can be written in the form 

M 



Pit) 



o 









c 10 c i |coi 

Vo o o 



cio| c 10 c 01 

2 







l-|cio| 2 -|c i| 



(12) 



The dynamics of the two qubits is therefore completely 
characterized by the two amplitudes: 



cio(i) 
coi it) 



r^\ 2 + \r^\ 2 £(t)]c w (0) 

r «M 2 )[l-£(i)]coi(0), (13) 
r (D r ( 2 )*[l-£(i)]c 10 (0) 

| r (i)|2 + | r (2)|2 f(t) l Coi(0)j (M) 



with r W = x U) /\xtI where \ X t\ = Vlx (1) I 2 + lx (2) I 2 , 
and 



£{t) = e -("-«*)*/4 



cos 



n g t 



K — l28 {Slqt 

■ sm 1 



2il r , 



(15) 



where tt g = -J 4\xt\ 2 W 2 + 5 2 + iSn — k 2 /4 is the gener- 
alized Rabi frequency. Note that Eqs. (|13[) and Q14[) have 



exactly the same structure as Eqs. (|2|) and (J3]), obtained 
for the single-mode Dicke model without losses. For- 
mally, the cavity losses appear as an additional imaginary 
part of the detuning 5 i— > 5 + ik/2. Accordingly, the effect 
of the cavity losses is described by the modification of the 
time-dependent coefficient £(t), which is now damped at 
rate k/4, and by the K-dependent shift of the Rabi fre- 
quency. For k — ► 0, the Lorentzian spectral density pop 
tends to Dirac's delta distribution, J{uj) W 2 8(w— wc), 
and Eq. (JTSJ reduces to Eq. (gj), with a ij) = x^'W. 

It is worth noticing that, as one sees directly from 
Eq. ©, the subradiant state 1^-), given by Eq. ([5]), is 
still decoupled from the vacuum cavity field. Hence, if the 
atomic system is initially prepared in this state, no ex- 
change of excitation with the cavity field will take place. 



III. EFFECTIVE MODEL OF ION-CAVITY 
INTERACTION 

A. Physical setup 

Ion-cavity QED experiments use calcium ions which 
are trapped in a linear Paul microtrap and interact with 
a quantized mode of a high-finesse optical cavity [1, . In 
Fig. [2] we show the relevant energy-level structure, cou- 
plings, and decay channels for the compound system of 
two 40 Ca + ions and a single cavity mode. The atomic 
ground state 4 2 S , 1 / 2 is coupled to the electronically ex- 
cited state 4 2 P!/ 2 by a (classical) pump laser injected 
from the side of the cavity. On the other hand, the ex- 
cited state 4 2 P!/ 2 is coupled to a metastable state 3 2 L> 3 / 2 
by the quantized cavity mode. The excited state 4 2 P 1 / 2 
decays spontaneously to the states 4 2 S , 1 / 2 and 3 2 D 3 / 2 at 
rates 7s and 7d, respectively, and the cavity photon is 
damped at rate explained in the previous section. 

A realistic theoretical description of the dynamics of 
a single 40 Ca + ion coupled to the cavity mode has been 
given in Ref. [3(|. The authors consider there also the 
effect of cavity losses and spontaneous emission, taking 
into account all the Zeeman sublevels of the three rel- 
evant electronic states. The main consequence of the 
presence of the Zeeman sublevels is a reduction in the 
coupling driven by the cavity field by a factor of \/3 with 
respect to the simpler three-level model considered here. 
Therefore, we will use in the following a three-level model 
scheme with such a reduced effective coupling to account 
for the presence of the Zeeman sublevels. In the experi- 
ments, the ions sit at the bottom of the trapping potential 
and arc cooled down to the Lamb-Dickc regime. Under 
these conditions one can assume that the ions are kept 
at fixed positions and neglect recoil during the emission- 
absorption process. 

In the following we will consider as initial atomic states 
those in which one of the two atoms is in its ground 
state and the other one is in its metastable state, i.e., 
the states \S^D^) and \D^S^). In order to prepare 
these staters, if the vibrational sidebands are not resolved, 
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FIG. 2: (Color online) The relevant electronic states of two identical Ca + ions and corresponding couplings provided by an 
external pump laser and a quantized cavity mode. The excited electronic state decays spontaneously to the ground and to the 
metastable states, and the cavity mode is damped as well. 



it is necessary to use a selective laser addressing of the 
individual ions. This is routinely done in trapped-ion 
experiments with 40 Ca + ions (see, e.g. 



respectively, 



311) 



The two identical ions interact with the quantized cav- 
ity mode of frequency ujc via laser-assisted two-photon 
processes, as shown in Fig. [5J The ions are irradiated 
by a laser beam of frequency lul + The laser beams 
and the cavity field are far detuned by A from the elec- 
tronic level |-P^), such that ujp — ujd = wc + A and 
ujp — ujs = ujl + A. Therefore, the setup provides each 
ion j = 1,2 with a Raman coupling between the levels 
\SV>) and \DW). 

The time evolution of the composite system of the two 
ions and the cavity mode can be described by a master 
equation 



= — i [if (t), £)] — ~ (aJ ag + ga) a — 2aga)} 



3=1,2 
7£ \ - 
2 ^ 

3=1,2 



A%g- 



gA<p P -2A$ P QA%), (16) 



where we have included the cavity field damping at rate 
k, the spontaneous emission channels (two for each ion) 
at rates 7s and jjy, and where the coherent dynamics is 
generated by a Hamiltonian 



H{t)=co c (a^a + 1 -)+Y J E ^ A U 

3 = 1,2 l=S,P,D 

E (e-^+s^tgU)^ + g^ a A% + H.c.) . 



3=1,2 



The atomic operators are defined as A\f) = 
with 1,1' = S,P,D and j = 1,2. Finally, the coherent 
couplings provided by the laser and the cavity mode are, 



(17) 



Z'Cj)| 



(18) 
(19) 



with fc/, and kc being the wave numbers of the laser and 
the standing cavity mode. 



B. Effective two-level model 



When the detuning A is sufficiently large compared to 
the couplings, A ^> , , the excited electronic states 
can be adiabatically eliminated from the dynam- 
ics, as described in detail in Appendix [XI In this case the 
system can be effectively described as composed of two 
two-level atoms interacting with a cavity mode. For this 
purpose, we denote the ground and the metastable states 
of the jth atom as = \S^) and |0 ( ^) = \D^) 

(N.B., the true atomic ground state corresponds to the 
excited state of the effective two-level system, since it is 
able to emit a cavity photon through the Raman transi- 
tion). 

The adiabatic elimination of the excited levels {l-P^-*)} 
is not at all trivial due to the inclusion of the spontaneous 
emission processes [32| . We show in Appendix [2] that an 
effective Tavis-Cummings Hamiltonian can be derived, 
describing an excitation exchange between the ions and 
the cavity. However, one needs to include (i) two Stark 
shift terms per ion (one, in particular, being dependent 
on the state of the cavity mode) and (ii) an overall re- 
scaling of both the free and the interaction energies by a 
factor explicitly dependent on the emission rates. 

It turns out that, in the interaction picture with re- 
spect to Hq — A J2j Ap P , where Hq is given by the first 
two terms on the right-hand side of Eq. (|17[) . the coherent 
part of the evolution of the ion-cavity system is described 
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by an effective Hamiltonian 



H, 



3=1,2 



A 



, A U) 
^11 



(20) 



where the position-dependent parameters ft^ are defined 



as 



ft) = e *k L x^ gin (j fcoa .(j)) j 



and the dimcnsionlcss rcnormalizing prcfactor is 



A ; 



A 2 + (7s + 7d) 2 /4' 



(21) 



(22) 



This Hamiltonian resembles the Tavis-Cummings Hamil- 
tonian (JTJ) , except for the photon-dependent Stark shift 
term. However, since the original microscopic model in- 
cludes dissipative processes, the unitary evolution gen- 
erated by -f/eff needs to be supplemented by decohering 
terms that have a very peculiar structure. Indeed, the ef- 
fective master equation that describes the time evolution 
of the ions and the cavity contains four (now both dissi- 
pative and non-dissipative) processes (described by jump 
operators) that take into account the effects of the spon- 
taneous emission as seen in the restricted atomic sub- 
spaces spanned by {|0^), The cavity damping 
appears in the restricted subspace in the same form as in 
the original model. 

The effective master equation reads 



dg 



— i [H e ff, g] — — (a^ag + go) a — 2aga^) 



3=1,2 
m=S,D 



(23) 



where the jump operators for each ion j are given by 

c ti) = e -iS L t Q A U) + ft)* g aA U) ; ( 24 ) 
e-^nA^+P^gaA^, (25) 



C 



D 



while the effective decay rates arc r„ = £,"/ m / A 2 , where 
m = S, D and the prefactor £ is given by Eq. (f2"2"| . The 
structure of these jump operators is easy to interpret once 
the full level configurations of Fig. [2] are taken into ac- 
count. Let us consider, for example, the operator C§ of 
Eq. ft arises from the spontaneous emission process 
^ 2 P\/2 — > 4 2 Si/2 of the jth atom, now being restricted to 
the two-dimensional subspace {|0w), The jump 

operator Co has two contributions, both of them de- 
scribing non-dissipative decoherence by pure dephasing 
processes (as one understands from the fact that they do 



not produce any excitation loss). These two contribu- 
tions account for the interruption of the ion-cavity exci- 
tation exchange (vacuum Rabi cycle) by the spontaneous 
emission. The first term is an unwanted repopulation of 
state occurring after the laser has virtually brought 
the system to the intermediate level \P^)) of the full Ra- 
man cycle. The second term is also due to decay into 
state but this time the virtual excitation of level 

|PW) is caused by the cavity field. In conclusion, both 
processes interrupt the vacuum Rabi cycle without the 
excitation being lost as, at the end, the two-level sys- 
tem is found in its excited state This implies that 
the excitation exchange can restart, but with a different 
phase. Thus, Cg describes a phase error. 

A similar interpretation scheme can be adopted for the 
two terms constituting C$ in Eq. (|25|) . However, this 
time the involved process is the spontaneous emission 



4' Pi 



1/2 



3 D%/2- Whether it occurs after the virtual ex- 
citation of level \P^) performed by the laser (first term) 
or by the cavity field (second term), the result is that at 
the end the two-level system is found in its ground state 
|0 ( ^) and that one excitation has been lost either from 
the atom or from the cavity mode. Therefore, this jump 
operator causes dissipative decoherence. We note that, 
at this stage, the four jump operators of Eqs. (|2"4")) - (|2"5|) 
are both explicitly time dependent and implicitly posi- 
tion dependent via the coefficients ft) [ see Eq. (|2T|)]. 

A phase rotation within the restricted Hilbert space, 
spanned by the states with at maximum one excitation, 
allows transforming the effective Hamiltonian (|20p into 
the Tavis-Cummings Hamiltonian |T]) as well as remov- 
ing simultaneously the time dependence from the jump 
operators (|24|) and (|25|) . This is described in Appendix IB"1 
Therefore, in a suitable rotating frame, the following ef- 
fective Tavis-Cummings Hamiltonian is obtained: 



rreff _ ; ,eff/_ f„ , 1 \ , , ,eff JJ) 

3=1,2 

+ E (oc^a^a^+a^a^), (26) 

3=1,2 

where we have introduced again the spin inversion opera- 
tors used in Sec.|TTJ The effective Dicke model parameters 
are 

ff _ 2|/3 Tff | 2 



Jc ~ s 3A ' 



(27) 

-ffr ( (W \PT9\ 2 \ ( . 

Ua-*l-Z{-£ (28) 



(29) 



where \(3 T \ = yWW+WW- Th e effective detuning 
is given by 

r , ,cff , ,off r t l^l 2 ~ \Pt9? 

o e s = uJ A -oj c = d L - £ , (30) 



7 



and the relative effective coupling strengths 

[cf. Eqs. and©] are directly given by the 

position-dependent parameters , since now = 

a^/\a eS , T \=f3 (j) /\f3 T \. 

Comparing Eqs. (J26]) and {(23j) with Eqs. (TJ and (fTTjl . 
respectively, we see that, when the effective atomic spon- 
taneous emissions are negligible, this system allows us to 
realize the Dickc model in the non-ideal cavity case. 



C. Effective spontaneous emission processes 

As mentioned before, we restrict our study to the case 
in which only one or zero quanta are present in the com- 
posite system of the two ions and the cavity mode. There- 
fore, the compound state of the two atoms and the cav- 
ity photon can be expressed in the basis (lO^O^O^), 
IfjWh^)), |0«l( 2 )0( c )>, JlWO^O^))} (see Ap- 
pcndix[B|. Consequently, the jump operators ([24j)- (J25J) 
can be normalized with respect to the operator norm 
| A\\ = sup||0|| =1 || A \<f>) || , where \<t>) belongs to the Hilbcrt 
space spanned by the basis defined above. The introduc- 
tion of the normalized jump operators allows to define 
the effective spontaneous emission decay rates Tm un- 
ambiguously. 

The normalized jump operators are 



ci x) = iiWo^xsii, 



4 2) = |o (1) i (2) o( c ))<$ 2 |, 
cg> = loWoVxsii, 

C f = loWoVX^al, 



where the decaying states are 



l*i) 
l* 2 ) 



n* iiWqWqc?)) + (3^g* |oWo( 2 )i( c )) 

ft* | (l)l(2) (^)> +/ 5(2) g * |0(l)0( 2 )l(O) 



VM 2 + \P {2) 9? 

The corresponding rescaled decay rates are given by 



{ S = i{M 2 



\0 u) 9\ 2 )f 2 , 

\^9\ 2 ff 2 - 



(31) 
(32) 
(33) 
(34) 

(35) 
(36) 

(37) 
(38) 



The cavity photon annihilation operator a = 
| (i)0(2) (cj^o( 1 )o( 2 h( c ')| is already normalized in 
our restricted basis. 

The spontaneous emission decay rates for the consid- 
ered states of a calcium atom are 7s/27r = 22.3 MHz 
and 7d/2tt = 1.7 MHz. Therefore, T { f > and the 

dominant effective spontaneous emission jump processes 

(7) 

are described by the operators C s : . Consequently, ac- 
cording to the discussion above, the main decoherence 
sources are the non-dissipative dephasing processes that 
conserve the energy of the ion-cavity system. 



The character of the decaying state, and hence the 
corresponding jump operator is defined by the balance 
between the strengths of the laser pumping O and the 
cavity coupling p^'g. In the strong laser pumping case 
(|fi| ^> \P^g\) the non-unitary dynamics of the atomic 
reduced system is dominated by phase diffusion processes 

described by the operators . In the weak laser pump- 
ing case (|f2| <C \P^'g\), on the contrary, the processes 
described by the operators aA^ dominate. Moreover, 
as one can see from Eqs. (|35|) and (|36"1) . one can further 
modify the character of the specific atomic jump oper- 
ators by changing the relative position of the ions with 
respect to the cavity field through the (3^ parameters. 

The significance of the spontaneous emissions can be 
estimated by the ratio 



r 



'cir 



1 + |/3(%/Q| 2 7s 
\pV)g/Sl\ A ' 



(39) 



For a fixed detuning A this ratio has its minimum value 
275/A when \/3^'g/fl\ = 1, i.e., when the couplings pro- 
vided by the laser and the cavity field are equally strong. 
On the other hand, for fixed coupling strengths, the ra- 
tio is inversely proportional to the detuning A. This can 
be exploited in order to minimize the role of the effec- 
tive spontaneous decay. The cavity damping k is neither 
affected by the detuning nor the couplings. 

Finally, we note that for large detunings, A ^> 75,71), 
the dimcnsionlcss prefactor £ ~ 1 and the effective de- 
cay rates as well as the effective coupling terms have 
simplified expressions. The effective couplings are then 
given by acQ ~ —j3^g*Cl/A, while in the limit of strong 
and weak laser pumpings the dominating decay rates are 



- |f2| 2 7 S /A 2 andr^ j - |/3(%| 2 7 S /A 2 , respectively. 



,U) 



IV. ENVIRONMENT-INDUCED 
ENTANGLEMENT: RESONANT REGIME 

In this section we study, analytically and numerically, 
the dynamics of the entanglement between the electronic 
degrees of freedom of the two atoms. The generation 
of entanglement between the ions and its persistence at 
long times are, indeed, a clear manifestation of the collec- 
tive (subradiant) behavior. In particular, entanglement 
generation is mediated by the interaction with the quan- 
tized cavity field which is initially prepared in the vac- 
uum state. If the atomic spontaneous emission processes 
are negligible and we face the bare Dicke model, the dy- 
namics can be described exactly. We compare these ex- 
act analytical results to numerical simulations including 
the spontaneous emission effects. The simulations were 
implemented by u sing the Monte Carlo wave function 
(MCWF) method f33T[34l| . We begin by considering the 
resonant case, where the effective detuning <5 ff = 0, with 
S oS given by Eq. J30]). 



8 



A. Analytical solution neglecting spontaneous 
emission 

The effective model describing the dynamics when 
spontaneous emissions are negligible is given by the mas- 
ter equation (jlip with the effective Tavis-Cummings 
Hamiltonian (|26|) . as described in Sec. Ill Bl The ana- 
lytical solution for the atomic density matrix is given by 
Eqs. (fT2|) and <fl5|l. with X U) W = a { £ = f3^g oS . 

We are interested in the collective dynamics when ini- 
tially one excitation is present in the atomic system and 
the cavity is in its vacuum state. Any initial atomic state 
containing one excitation can be written in terms of the 
supcrradiant and subradiant states ©-{SI) as 

Mo)> = /? + |v+>+/3-hM- ( 4 °) 

As time passes, the collective atomic state decays via the 
evolution of the superradiant component, 

w+im) =£(t)p+, (4i) 

with £(t) given by Eq. (|15[) . The subradiant component 
(ip-\ip(t)) = f3—, however, remains unchanged. Conse- 
quently, for times, such that nt ^> 1, the atomic state 
will be in general a statistical mixture of the collective 
ground state lO^W 2 )) and the subradiant state with 
weights dependent on /3__ , which in turn depends on the 
relative coupling strengths rd> . 

In the following we focus on the dynamics of entangle- 
ment between the atoms. In order to quantify the sta- 
tionary asymptotic entanglement of the final state we use 
Wootters's concurrence [35| which, for a density matrix 
of the form of Eq. (112[) . is given by 

C(t)=2\c w (t)c* 01 (t)\, (42) 

with cio(t) and c i(t) given by Eqs. (fT^|) and (TT4")) . In 
general, the concurrence is zero for factorized states and 
unity for maximally entangled states. For > 1 we 
obtain a stationary concurrence value 

C stat = 2|r( 1 V( 2 '||/3_| 2 . (43) 

As expected, the value of the stationary concurrence is 
directly related to the subradiant component of the ini- 
tial state. If both atoms arc coupled to the EM field, 
the stationary value of the concurrence, for any initial 
state with /3_ ^ 0, will be nonzero. When the atoms are 
initially prepared in the superradiant state, i.e., /3_ = 0, 
the system approaches asymptotically the pure factorized 
state |0«0( 2 )). 

For the initially factorized states [l^O^ 2 - 1 ) and 
|0 (1) 1< 2 )), the interaction with the environment gener- 
ates entanglement in the atomic system. For these initial 
states the stationary concurrence takes the values C s t a t = 
2| r «|(i - | r (i)|2)3/a and Cstat = 2|r«£Vl-|r(i)| 2 , re- 
spectively. As we have noticed in Ref. [29[ , the factorized 
states are those that maximize the stationary concur- 
rence for certain values of r' 1 '. The maximum value of 



stationary concurrence, for both the two factorized initial 
states considered here, is C^ff = max i r (i)i g ro .1] Cstat — 
0.65. This value is obtained with \r^'\ = 0.5 and 
|rW| ~ 0.87 (i.e., \r^ \ = 0.5) for initial states jl^O^) 
and lO^l^ 2 )), respectively. 

We note in passing that when only one of the two atoms 
is coupled to the EM field, i.e., r-W = or r' 2 ' = 0, 
the stationary concurrence is zero. In this case, indeed, 
the subradiant and the supcrradiant states coincide with 
states Jl^O^) and |0^ 1 - ) 1^ 2 ^ as one can see from defi- 
nitions (|5|) and ([6]). 

From the definition of the generalized Rabi frequency 
given by Eq. (|15p . which in the resonant case reads as 
il g = ^/4|/3x5 G ff| 2 — k 2 /4, two extreme regimes can be 
defined. In the weak ion-cavity coupling regime, defined 
by 4|/3<r(7 ff I ^ k j the generalized Rabi frequency is purely 
imaginary. Therefore, according to Eq. (p~5j) . the Dickc 
model predicts a solution given by monotonic hyperbolic 
sine and cosine functions. The opposite limit is the strong 
ion-cavity coupling regime, defined by 4|/3T5cff | 3> k- In 
this case the generalized Rabi frequency is real and the 
Dickc model predicts damped oscillatory dynamics. 



B. MCWF simulations in the presence of 
spontaneous emission 

In this section, we focus on the effect of the spon- 
taneous emissions on the subradiant-state-based en- 
tanglement generation described in the previous sec- 
tion. We consider again as initial atomic state 
\tp(0)) = \1^0^) with the cavity in the vacuum 
state |0( c )). For a given value of r^ 1 ' G [0,1], 
we choose (3^ and (3^ to be positive real numbers 
such that the larger of the two is always unity and 
the smaller one is minjr^ 1 ^ /Vl — rW 2 , Vl — r( 1 ) 2 /r' 1 '} 
[cf. definition Now \(3 T \ 2 = |/? (1) | 2 + |/3^| 2 = 

min{l/rW 2 ,l/(l - r^ 2 )} € [1,2]. The physical param- 
eters have been chosen in accordance to the experiments 
of Ref. H and are summarized in Table |U The size of the 
ensemble in the MCWF simulations is N — 1000. We are 
using the variant of MCWF method described in [34| . 

The value of the cavity coupling constant g in Table Q] 
refers to the new miniature trap recently realized at the 
University of Sussex [36[ . The reference value kq for the 
cavity damping can nowadays be improved by at least 
one order of magnitude. Finally, the detuning A can be 
easily increased in the experiments, with respect to the 
reference value Ao. 

With the experimental parameters of Tabic Q] the cou- 
pling strengths O and g are of the same order. Therefore, 
neither the strong nor the weak laser pumping regimes, 
introduced in Sec. IIII C[ are reached and, consequently, 
all the effective decay processes caused by the sponta- 
neous emission arc combinations of two different physical 
operations, as interpreted in Sec. IIII Bl 

Let us denote the atomic density-matrix components 
as p abtCd ee (aW&^HcW^ 2 )}, where a,b,c,d = 0,1. 
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TABLE I: Values of physical quantities used in the simula- 
tions. Note that the cavity coupling is here explicitly scaled 
by the Clebsch-Gordan factor l/\/3 and, in the text, also by 
the position-dependent parameters /3^'. 



Quantity Symbol Value (2tt MHz) 



Laser coupling 


Q 


9.0 


Cavity coupling 


9 


6.5 /VS 


Decay rate 4 2 P 1/2 -> 4 2 S 1/2 


is 


22.3 


Decay rate 4 2 P 1/2 — > 3 2 D 3 / 2 


JD 


1.7 


Detuning 


A 


20.0 


Cavity damping 


no 


1.2 



the effective spontaneous emissions in favor of the coher- 
ent dynamics, as explained in Sec. If 11 CI The situation 
is clarified in Fig. [3] which shows the scaling of the effec- 
tive coupling strength g c g and the dominant spontaneous 

emission decay rate [cf. Eqs. (J29]), ([37]), and ([38])] as 
functions of detuning A. The cavity damping rate k is 
not affected by the detuning. The relative values of the 
three key parameters g e g 7 Tg , and k characterize the 
dynamical regime: (i) the ratio \g e g\/n defines the strong 
and the weak ion-cavity coupling regimes and (ii) the 
magnitude of compared to \g e g\ and K, in turn, de- 
scribes the significance of the spontaneous emission pro- 
cesses and tells us whether the dynamics is well described 
by the Dicke model or not. 



N 

x 



> 

o 

z 
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o 

LU 
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10 
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FIG. 3: (Color online) Scaling of the effective ion-cavity cou- 
pling \g c s\ (middle line for large A) and effective sponta- 
neous emission decay rates Fg (lowest line for large A; with 
= 1) as a function of the detuning A. The isocurves 
4\Prg c g\/ k. — const (thin lines) are parallel to the \g e g\ curve, 
so that the weak ion-cavity coupling regime is in the upper 
right corner and the strong ion-cavity coupling regime in the 
lower left one. The cavity decay rate k (horizontal line) does 
not depend on the detuning. The effective spontaneous emis- 
sion events are suppressed for large detunings. 



The density matrix remains still in the same block form 
of Eq. (II2|) even in the presence of spontaneous emissions. 
The concurrence is therefore given by C(t) = 2|/9oi,io(i)|- 
In the following we will examine the effect of the 
spontaneous emissions by comparing the concurrence 
as a function of time for fixed values of 4\PT9es\/ K = 
4|/?t£s , ^/ k A|. We study large detunings (A 3> Is, Id), 
so the prefactor £ ~ 1. In the examples we change k and 
A, such that k/k = 0.1, 0.01 and A/A = 10, 100, 1000, 
while keeping the product kA constant. Physically, this 
corresponds to using different cavity qualities and de- 
tunings which, furthermore, influences the effective dy- 
namical parameters. Larger detunings, indeed, suppress 



1. Weak ion-cavity coupling regime 

In this regime, the oscillatory dynamics stemming from 
the coherent coupling between the atoms and the cavity 
is heavily damped. In Fig. 3] we plot the concurrence 
as a function of both time and the relative coupling 
strength n 1 ' for A/Arj = 100 and k/ko = 0.1, giving 
| <7cff | / 27r = £<7f2/27rA = 17 kHz. All the other param- 
eters are chosen as in Table [TJ We recall that initially 
the atomic state \tp(0)) = |1 (1) (2) ) is factorized. The 
initial dynamics of the concurrence shows a monotonic 
increase, as the supcrradiant component [sec Eq. (|40[) ] 
rapidly fades away while the subradiant component re- 
mains intact. However, because of the presence of sponta- 
neous emission, the subradiant state is not anymore per- 
fectly decoupled from the dynamics and, consequently, 
the concurrence will not reach a steady-state value. 

We note that the best peak value of the concurrence 
C ~ 0.6 is achieved for ~ 0.55, i.e., as expected, for 
an asymmetric configuration (r^ 1 ' ^ 1/V%) of the ions 
with respect to the cavity field. However, this value of 
n 1 ' is now slightly different than the one obtained in 
Sec. I IV Al where spontaneous emissions were neglected 
(jA 1 ) = 0.5). We will further discuss this point when con- 
sidering the position dependence of the jumps statistics 
at the end of this subsection. 

In Fig. [5] we further study the effect of the spontaneous 
emissions in the weak ion-cavity coupling case. In this 
figure, we compare the predictions of the Dicke model, 
described in Sec. Ill Bl with the dynamics of the ion-cavity 
system in the presence of the spontaneous emissions for 
A/A = 100, k/k = 0.1 and A/A = 1000, k/k = 0.01. 
The dynamics of the concurrence clearly shows that, in 
the first case (k/kq = 0.1), the system approximates the 
Dicke model well while \fl g \t/2ir < 2.5, where the gen- 
eralized Rabi frequency \Sl g \ is given by Eq. (fT5]) . For a 
better cavity (k/kq = 0.01), the concurrence approaches 
its quasi-stationary value and the system approximates 
the ideal Dicke dynamics for longer times, |f2 s |t/27r < 20. 

We finally look at the statistics of the quantum jumps, 
described by the jump operators Cm in Eqs. ([3"Tj) - (|3"4"]) . 
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FIG. 4: (Color online) Concurrence as a function of time 
and the relative coupling strength r^ 1 ' in the weak ion-cavity 
coupling regime. The dynamics is initially monotonic since 
the existing superradiant component decays rapidly com- 
pared to other dynamical time scales. The subradiant state 
component decays eventually because of the atomic spon- 
taneous emissions. The best entanglement production oc- 
curs with asymmetric couplings (r' 1 ' 7^ l/\/2). Parameters: 
A/Aq = 100, k/k = 0.1 



FIG. 5: (Color online) Dynamics of the concurrence in the 
weak ion-cavity coupling regime for r (1) = 0.55. In the Dicke 
model with cavity losses (highest line) a stationary value of 
the concurrence is reached as the superradiant component is 
over-damped. Parameters: A/Ao = 100, k/ko = 0.1 (lowest 
line, 2iv/\fl g \ = 23 /is); A/A = 1000, k/k = 0.01 (middle 
line, 2n/\fl g \ = 230 fis). 



First of all, wc note that the source states [sec 

Eqs. ((33) and $o§] of the jump operators and C$ 
are identical for a given atom j = 1,2. Therefore, the 
jump statistics of the two corresponding decay channels 
will also be the same with a branching ratio given by 

r^/To = Is/lD ^ 13. Our MCWF simulations con- 
firm that the dominant jump processes are those corre- 
sponding to the effective spontaneous emission operators 

(7) 

Cg and the cavity photon annihilation operator a. In 
Fig. [5] wc plot the average cumulative number of quan- 
tum jumps per ensemble member for the jump operators 

c s ] > c f\ and o. 

Looking at the statistics helps us to understand how 
the reservoir-mediated entanglement generation process 
depends on r' 1 ). We notice that the jump statistics of 
processes originating from the spontaneous emissions of 
atoms 1 and 2 are different. This is of course due to 
the asymmetry in the initial condition. Since initially 
the excitation is present in atom 1, the average cumula- 
tive number of jumps per ensemble member is typically 
greater for Cg than for . The peak in the cumula- 
tive number of jumps, for the three different jump oper- 
ators considered in Fig. [5J moreover, is reached in corre- 
spondence of different values of r' x >. This indicates that 
the value ~ 0.55, which optimizes the concurrence 
generation (see Fig. [4J , corresponds to a compromise be- 
tween the different A 1 ' -dependent jump statistics. In 
particular, the deviation from the optimal value in the 
absence of spontaneous emission (r^ = 0.5) might be 




TIME his) 



FIG. 6: (Color online) Average cumulative number of quan- 
tum jumps per ensemble member for each decay channel in 
the weak ion-cavity coupling regime. From above: C s 1 , C s 2 \ 
and a. Parameters: A/Ao = 100,k/ko = 0.1. 



due to the fact that the number of C s -jumps increases 
for decreasing values of 
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FIG. 7: (Color online) Concurrence as a function of time 
for different values of the relative coupling strength r' 1 ' in 
the strong ion-cavity coupling regime. Oscillations appear 
because of a relative phase evolution between the superradiant 
and subradiant states. Parameters: A/Ao = 10 and k/ko = 
0.1. 



2. Strong ion-cavity coupling regime 

In the strong ion-cavity coupling regime, the cavity 
damping is slow compared to the coherent dynamics. 
Therefore, a slowly damped oscillatory behavior of the 
concurrence is expected. In Fig. [7] we plot the concur- 
rence as a function of both time and the relative cou- 
pling strength A 1 ' for A/A = 10 and k/ko = 0.1, giving 
\g c g\/2ir — £<?r2/27rA = 170 kHz. All the other parame- 
ters are chosen as in Tablefl] Note that the ratio \g c g\/K 
is now one order of magnitude bigger than in Sec. lIVBTI 
The dynamics has an oscillatory character, since the su- 
perradiant component survives much longer than in the 
weak ion-cavity coupling regime. However, due to the 
presence of the spontaneous emissions the concurrence 
does not reach a steady-state value in this regime either. 

The best peak value of the concurrence, C ~ 0.6, 
is now obtained for A 1 ^ ~ 0.46. In Fig. [5] wc choose 
this value of A 1 ' and we compare the dynamics of the 
single-mode Dicke model with cavity losses to the dy- 
namics of the ion-cavity system in the presence of effec- 
tive spontaneous emissions for the cases of A/Ao = 10 
with k/k = 0.1, and A/A = 100 with k/k = 0.01. In 
the second case, i.e., for a better quality factor, the sys- 
tem approximates the Dicke model for longer time scales, 
as one would expect. In this case one can clearly observe 
the damped Rabi oscillation at the generalized Rabi fre- 
quency given by Eq. ((15)1. 

It is worth noticing that, in the strong ion-cavity cou- 
pling regime, the laser-mediated interaction with the cav- 
ity vacuum allows us to generate a highly entangled state 
of the two ions, as one can see in Fig. [5J In particular, 
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FIG. 8: (Color online) Time evolution of the concurrence 
in the strong ion-cavity coupling regime with relative cou- 
pling strength r"' = 0.46. For the Dicke model with cav- 
ity losses (highest line), the concurrence approaches a con- 
stant stationary value after strong oscillations caused by 
the slowly decaying superradiant component. Parameters: 
A/Ao = 10, k/k = 0.1 (lowest line, 27r/|n 9 | = 2.7 /is); 
A/A = 100, k/ko = 0.01 (middle line, 2Tv/\n g \ = 27 fjs). 

for A/Ao = 100 with k/kq = 0.01, using a laser pulse 
of duration t ~ 27r/|fi g |, the generated state is close to a 
maximally entangled Bell state. 



V. ENVIRONMENT-INDUCED 
ENTANGLEMENT: DISPERSIVE REGIME 

In the previous section we have seen that by placing 
the ions properly, i.e., by adjusting the relative coupling 
strength A x \ it is possible to optimize the reservoir- 
mediated entanglement generation. The examples dis- 
cussed above deal with the resonant effective model, 
which is defined by the condition 8 c s = 0, which in 
turn corresponds to a physical detuning 6l = £[^ 2 — 
(/?ts0 2 ]/A [cf. Eq. ([30])]. We have seen that the highest 
value of the concurrence is obtained in the strong ion- 
cavity coupling regime. 

In Ref. (37[, however, the single-mode Dicke model 
with cavity losses is studied in the dispersive regime, 
showing that a high degree of entanglement can be ob- 
tained also in the weak ion-cavity coupling regime. For 
this reason we now look at the off-resonant entanglement 
generation process in the ion-cavity QED, i.e., we con- 
sider the case in which <5 G ff ^ 0. In the dispersive regime, 
the relative position of the ions does not play an essen- 
tial role and in fact one shows that the optimal value of 
A 1 ' is obtained for equal coupling of the two ions, i.e., 
r«=r( 2 ) = 1/V2 M- 

We consider once more the initial atomic state IV'(O)) = 
|l(!)o( 2 )) combined with the cavity in the vacuum state 
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FIG. 9: (Color online) Concurrence as a function of time for 
different values of the detuning Sl for a homogeneously cou- 
pled case (r*- 1 ' = r^ 2 ' = l/s/2). The entanglement generation 
slows down when passing into the dispersive regime. The 
resonance is dislocated from the origin because of the Stark 
shifts. Parameters: A/Aq = 10 and k/kq = 0.1. 

|0 (c) ). Wc set r-W = r (2) = l/y/2 (by choosing max- 
imally strong cavity-driven couplings = (3^ = 1), 
A/Ao = 10, and k/kq = 0.1, corresponding to the weak 
ion-cavity coupling regime of Sec. lIVBTI We now look at 
the time evolution of the concurrence for different values 
of the laser detuning Sl- Figure [9] shows the concurrence 
as a function of both time and detuning Sl ■ One can see 
clearly that the Stark shift terms appearing in the effec- 
tive Hamiltonian of Eq. (j2"0)) relocate the resonance con- 
dition from the origin to 5 l /2tt = ([Q 2 - (/3 T g) 2 ]/27rA = 
120 kHz. Figure [5] also shows that selecting the detun- 
ing Sl further away from the resonance produces higher 
values of concurrence. In particular, with the chosen pa- 
rameters the maximum value of concurrence C ~ 0.62 is 
obtained with Sl/2tt ~ 600 kHz. 

As demonstrated in Ref. [37| , increasing the detuning 
| (J e ff | correspondingly increases the time it takes for the 
concurrence to reach its peak value. The longer is the en- 
tanglement generation time, however, the stronger is the 
effect of the spontaneous emissions. In other words, the 
achieved gain in the entanglement generation obtained 
by increasing the effective detuning is quickly suppressed 
due to the spontaneous decay, as the overall time of the 
entanglement generation process increases. The maxi- 
mum value of entanglement achievable in the dispersive 
regime is therefore determined by the interplay between 
these two effects. 

It is worth noticing that going from the resonant into 
the dispersive regime changes the character of the gen- 
erated entangled state as well. To illustrate this point, 
we plot in Figs. \W\ and [IT] the populations and coher- 
ences, respectively, of the reduced atomic density matrix 
versus time and detuning Sl- These plots confirm the 




5 10 15 

TIME (lis) 



FIG. 10: (Color online) Populations of the atomic states 
poo, oo, Poi.oi, and pio,io (from above) as a function of time for 
different values of detuning Sl- Parameters are as in Fig. [9] 

increase in the entanglement generation time when going 
deeper and deeper into the dispersive regime (|<5 e ff| > 0). 
If we then focus on the dynamics of the coherences and, 
in particular, on the real and imaginary parts of the only 
nonzero off-diagonal element poi.io, we see that on reso- 
nance the imaginary part vanishes in accordance with the 
predictions of Sec. IIV1 Therefore, in the resonant regime 
the generated entangled state approximates the subradi- 
ant state. On the other hand, in the dispersive regime 
Rc[poi,io] — and Im[poi,io] 7^ 0. Indeed, in the absence 
of the spontaneous emissions, the generated state in the 
dispersive regime would be (|lWo( 2 )) ± i^l^)) /V2 
(positive sign for negative 5 e ff and vice versa). 

VI. SUMMARY AND CONCLUSIONS 

In this paper we have investigated how the single- 
mode Dicke model can be realized under experimentally 
feasible conditions using two trapped 40 Ca + ions inside 
a high-finesse optical cavity. We have taken into ac- 
count the spontaneous emissions of the ions as well as 
the damping of the electromagnetic field inside the cav- 
ity. In particular, we have derived an effective two-level 
description of the three-level ions interacting with the 
cavity mode. 

We have shown that under suitable conditions the two 
ions indeed behave collectively, with a coherent dynam- 
ical evolution well described by the Dicke model: two 
effective two-level systems exchanging an excitation with 
an effective one-dimensional cavity mode. The presence 
of decohering processes, such as the atomic spontaneous 
emission or the cavity field damping, modifies this ideal 



13 




5 10 15 



TIME (us) 



FIG. 11: (Color online) Dynamics of Re[poi,io], Im[poi,io], 
and |poi,io| (from above) as a function of time for different 
values of detuning Sl (observed concurrence is given by C — 
2|poi,io | ) - Parameters are as in Fig. 



nolo gy u sed in the context of ion-cavity QED experiments 
[1. 19U301] . In order to detect the generated entanglement, 
the state tomography of the atomic systems is needed. In 
recent years, this has been routinely performed in simi- 
lar trapped-ion systems, e.g., in the context of quantum 
computation and measuring the quality of quantum gates 
[3lj . Therefore, we expect our proposal to be within the 
reach of the experimental community. 
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picture. However, in the effective model, the spontaneous 
emission decay rates are proportional to 1/A 2 whereas 
the ion-cavity couplings scale as 1/A, where A is the 
detuning of the physical cavity frequency from the elec- 
tronic transition that it is driving. This difference in the 
scaling can be exploited in order to partly suppress the 
destructive effect of the atomic spontaneous emissions. 

We have identified the generation of entanglement as a 
fingerprint of the cooperative atomic behavior and ana- 
lyzed this process in detail. In particular, we have proven 
that it is possible to enhance the entanglement generation 
process by positioning the ions appropriately at different 
locations with respect to the standing mode of the elec- 
tromagnetic field inside the cavity. In the resonant case, 
where the two-level systems and the cavity mode have 
the same frequency, we have shown that asymmetric cou- 
pling with the cavity mode produces the highest degree 
of entanglement, even in presence of spontaneous emis- 
sions. We have studied both the weak and the strong 
ion-cavity coupling regimes, defined by the strength of 
the ion-cavity excitation exchange compared to the cav- 
ity field damping rate, and found out the optimal condi- 
tions for entanglement generation in both cases. 

Another possibility to optimize the entanglement gen- 
eration is to go to the dispersive regime in the ion-cavity 
coupling by using an off-resonant Raman transition. The 
maximum degree of entanglement in the dispersive and 
in the resonant regimes, for realistic values of the param- 
eters, is similar. Our results indicate, however, that the 
character of the generated entangled state in the disper- 
sive regime changes compared to the resonant case. 

Our experimental proposal is based on existing tech- 



APPENDIX A 

In this appendix we show how the master equation (JT5J) 
for A-coupled three-level atoms and a cavity photon 
is transformed into an effective two-level master equa- 
tion (|23[) by adiabatic elimination of the excited atomic 
states {\P^')}- Especially, the elimination transforms 
the jump operators related to the spontaneous emissions 
into the form given by Eqs. (j2"4")l and ([23]) . For the sake of 
generality, our treatment here is valid for N atoms and 
we allow each ion j to be irradiated by a separate pump 
laser with frequency u>l + $l and coupling strength gr/ . 

Passing into interaction picture g 1— » g = e lKt ge~ 
with respect to 

K = u c (<Ja+^+Y J £ <*i A u- a1 E A pp ( Al ) 

j l=S,P,D j 

transforms the operators as 

a i-> 5 = e -' lUct a, (A2) 

A%^Af D = e +i »° t Af D , (A4) 

A % - A^s = (A5) 
Af~Af=A$, l = S,P,D, (A6) 

for all atoms j, and the Hamiltonian (jTTJ) becomes 

j j 

(A7) 
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The dissipator part of the master equation (|16|) is invari- 
ant under this transformation. 

Let us define a projection V to a subspacc spanned by 
the to-be-eliminated atomic states {|P^)}, and another 
projection Q to the complementary subspace by 



0) 
pp 



Q 



(A8) 
(A9) 



Setting £iqp = qpq = 0, assuming qqq > £ PP , and 
neglecting the cavity damping in Eq. (|A15[) gives an ap- 
proximation 



Correspondingly, the density matrix divides into four sec- 
tions 



§ = QqQ + QgV + V~qQ + VgV 

= QQQ + QQP + QPQ + QPP- 

The Hamiltonian is similarly divided in parts 
Hpp = VHV = A pp, 



(A10) 



(All) 



QQP 



Qpq 



-A + z( 7s +7-d)/2 
A 2 + (7s+7,d) 2 /4 



QQqHqp- (A17) 



LQP = VHQ 



Hpq — H, 

= 52{gU>e-*S'*AU> 8 + g $aA 

3 

H QQ = QHQ = 



PD I ' 



(A12) 



(A13) 



We proceed to deriving an effective master equation 
for qqq, which describes the dynamics of a collection of 
effective two-level atoms and a cavity mode. Applying 
Eqs. (|A10j) - (|A13[) to master equation (JTHJ) gives 



QQQ 



iHqpQpq + iQqpHpQ - - ({a f a, Qqq} 



2aQ QQ ai) + 7s a sIqpp a 



PS 



Id 



J2 A dpQpp a 



PD' 



(A14) 



QQP - Qpq 



IS +7-D 

2 



- iA qqp + i 



QP 



iHqpqpp - - y{a ] a, qqp} - 2aQ QP a t J , 



(A15) 



Qpp =-{jS+ 1d)qpp - iHpqqqp + iqpqHqp 



- [{a^a, Q P p} - 2aQ PP a) 



(A16) 



Similarly, setting qpp = in Eq. (|A16[) and using the 
above approximations for qqp and qpq gives 



Qpp 



7s + 7-D 
1 



{HpqQqp - QpqHqp) 



Finally, by inserting Eqs. (|A17|) and (|A18|) into 
Eq. (|A14|) . we arrive at an approximated master equa- 
tion 
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QQQ 



-i 



(is + 7r>) 2 / 4 

7s 

^2A2 + ( 7S+7D )V4 



HqpHpq, Qqq 



1 {a f a, qqq} - 2aQ QQ a ] 



iO") < 

\PS* 



7d 



2A 2 +(7s+7d) 2 /4 



{(FQpA ps )(yl^],ffpQ), Qqq) - 2(A^ p Hpq)Qqq(Hq P A 

{(Hq P A%)(A%H pq ), Qqq} - 2(A% ) p H P q)bq Q (H qp A%) 



k r U} r 



j m=S,D 

r 



, (A19) 



which has the form of the master equation ([23]) . We can 
now recognize the effective Hamiltonian [cf. Eq. ([20)) ] as 



H, 



eft 



A 

A 2 , ( 7S+7 D -p H QP H PQ 



|„0')|2 I „(•?') 1 2 



(A20) 



the effective spontaneous emission jump operators 
[cf. Eqs. (ED) and J25])] as 



eg =AV p H PQ =g^e-^A^+g^aA%, (A21) 



Aj)-is[ j h A U) 



c 



4&#pq = S^e^'tAW+^oAW, (A22) 



(N.B., the operators are unique up to a global phase fac- 
tor e %9 ™\ where 9 m G K), and the corresponding decay 
rates as 



7s 



-.0') 
p 



A 2 + (7s + 7p) 2 /4 

TP 

A 2 + ( 7S + 7 p) 2 /4 



(A23) 



i 



To 
A 2 ' 



(A24) 



In the above equations the dimcnsionless prefactor £ is 
as defined in Eq. ([22]). 



APPENDIX B 

In this appendix, we review how the effective two-level 
Hamiltonian (|2"D|) with Stark shifts is matched exactly 
with the Tavis-Cummings Hamiltonian |T]) by passing 
into a rotating frame. Moreover, we show how to ex- 
ploit the same phase transformation in order to simplify 
the propagator for the numerical simulations. As in Ap- 
pendix [X] we assume each atom j = 1, . . . , N to be ad- 
dressed by its own pump laser with independent detun- 
ings 8^) and couplings g^ . 

We will restrict ourselves to the subspace with at most 
one excitation. We introduce the following notation for 



the atomic states: 

|^) = |0«.--0W), (Bl) 
1^')) = loW.-.l^-.-oW), (B2) 

so that the Hamiltonian ([20)) reads 

H =E [(e-^'A^bo 1 {C] )U 3) (c) | + ff.c.) 
j 

+ ^)|v ? ol (C) >(^ol (C) |, (B3) 

with Stark shifts = -^T, j \9c ) \ 2 / A and sU) = 

— £\gjj?'\ 2 /A, and effective ion-cavity couplings \W> = 

~^9c iSl I ' A - The time dependence of the coupling 
terms will be eliminated by a phase rotation of the basis 
vectors 



I^o0^> 
lv?ol (C) ) 



0') (c)) 



(B4) 

(C) ), (B5) 
qM), (B6) 



where /i, ^ G M are up to now free parameters. The 
Hamiltonian transforms accordingly into 

H^H'=J2 \<Po 1 (C) > 0< c > | + ff.c.) 

+ (5( c )+^)|^ol (C) )(^ol (C) | 

+ A ^o0^)(^o0( c )| (B7) 



The requirement of a full compatibility with the Tavis- 
Cummings Hamiltonian ([1]) within our restricted Hilbcrt 
space demands that /j, = /j,(i>) = (S^ + v)/3. Conse- 
quently, the effective Dickc model parameters arc identi- 
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ficd as [cf. Eqs. ([27j)-(p9]l] 

^ ) =^" ) +5W-^) + ^, (B9) 
ag=A^. (BIO) 
Moreover, the detunings are [cf. Eq. (|30l) ] 

«5« = _ w eff = S V) + S (3) _ (BU) 

This means that the resonance condition of the Dickc 
model (cf. Sec.lIVp is achieved with laser detunings 6^ = 
5(C) _ ^ wm i e i n th e dispersive regime (cf. Sec. [Vj) 

5^ 7^ 5^' — As a conclusion, having each (iden- 

tical) atom driven by their own pump laser allows us 
to simulate the inhomogeneous Dicke model ([!]) with in- 
dependent two- level transition frequencies lo^ an d cou- 
plings a»' . On the other hand, if one has only a single 
laser driving all of the atoms, the transition frequencies 
are the same for every atom j, but the coupling con- 
stants remain independent because of the position- 
er) 

dependent cavity couplings g£ . 

Another aspect of the performed phase transforma- 
tion is provided by how they affect the dissipator part of 



the master equation (|23|) . Within our restricted Hilbcrt 
space, the jump operators transform now as 

Cf 1 — > e-^+^'l^ lc) ){*j\, (B12) 
C U) ^e-^+^l^o^X*,-!, (B13) 
ai-^e-^M-^l^oO^)^!^! (B14) 

(global phase factors can be discarded immediately), 
where the decaying (un-normalized) states are 

\*i) = e^g% > \<p¥ ) 0W)+g%>*\<p l< c) )- (B15) 

In the numerical MCWF simulations, the dynamics is 
generated by a non-Hermitian Monte Carlo Hamiltonian 
H M c = H - |5Zm AmJlnJm, where H is the Hermi- 
tian Hamiltonian of the master equation, and J m and 
A m are all the jump operators and corresponding decay 
rates picked up from the dissipator part of the master 
equation. From the practical point of view, it is advan- 
tageous to have a time- independent Hmc-, since then the 
Dyson series of the propagator simplifies into exponen- 
tial form U(t,to) = exp[—iHMc(t ~ to)]- This is now 
achieved simply by choosing v = 0, and hence the phase 
transformation is unique. 
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